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Enhancements of Andreev conductance induced by the photon/vibron scattering 
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We analyze the subgap spectrum and transport properties of the quantum dot embedded between 
one superconducting and another metallic reservoirs and additionally coupled to an external boson 
mode. Emission/absorption of the bosonic quanta induces a series of the subgap Andreev states, 
that eventually interfere with each other. We discuss their signatures in the differential conductance 
both, for the linear and nonlinear regimes. 
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I. INTRODUCTION 

The bosonic modes, like photons 1.] or vibrational de¬ 
grees of freedom , can strongly affect electron tunneling 
through the nanoscopic systems [§]. When a level spacing 
of nanoobject is large in comparison to the boson energy 
wo and a line-broadening is sufficiently narrow, a series of 
the side-peaks Q may appear due to emission/absorption 
of the bosonic quanta. Such features (spaced by wo) have 
been really observed in measurements of the differential 
conductance for several nano junctions [BUI- 

Similar bosonic modes are currently studied also in the 
systems, where the quantum dots/impurities are coupled 
with superconducting reservoirs [9H20[. Since the prox¬ 
imity effect spreads electron pairing onto these quantum 
dots, the bosonic features manifest themselves in a quite 
peculiar way. They could be observed by the Josephson 
@4111 and the Andreev spectroscopies |12h17| , in photon- 
assisted subgap tunneling [l8j, transient phenomena 0, 
or in prototypes of the nano-refrigerators operating due 
to the multi-phonon Andreev scattering |20j |. 

First of all, in a subgap regime (assuming wo smaller 
than energy gap A of superconductor) the bosonic fea¬ 
tures are expected to be more numerous than in the nor¬ 
mal state. This is a consequence the proximity effect, 
mixing the particle and hole excitations. Secondly, it has 
been shown numerically flB . flB . |20| that the linear (zero- 
bias) Andreev conductance exhibits the bosonic features 
spaced by a half of wo- To our knowledge, this intrigu¬ 
ing theoretical result was neither clarified on physical ar¬ 
guments nor checked experimentally. Verification would 
be feasible by the tunneling spectroscopy using e.g. low- 
frequency vibrations of some heavy molecules or slowly- 
varying ac electromagnetic field. Let us emphasize, that 
such low-energy boson mode need not be related with 
any pairing mechanism of the superconducting reservoir. 

The purpose of our paper is to provide a simple ana¬ 
lytical argument, explaining the reduced frequency w o/2 
of the bosonic features in the linear Andreev conduc¬ 
tance versus the gate-voltage. We also study in detail the 
multiple subgap states originating from the boson emis¬ 
sion/absorption processes. We analyze their signatures 
both in the quantum dot spectrum and the tunneling 
transmission. The latter quantity can be probed by the 
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FIG. 1: (color online) A scheme of the quantum dot between 
the metallic (N) and superconducting (S) electrodes and cou¬ 
pled to the monochromatic boson (phonon or photon) mode. 


(low-temperature) differential conductance as a function 
of the source-drain bias. We predict that the multiple 
Andreev states could be seen with a period, dependent 
on the gate voltage. 


For calculations we consider the setup displayed in fig¬ 
ure [I] It can be practically realized in a single elec¬ 
tron transistor (SET) using e.g. the carbon nanotube 
suspended between the external electrodes (like in Refs 
[5, 6|). Another possibility could be the scanning tun¬ 
neling microscope (STM), where the conducting tip (N) 
probes some vibrating quantum impurity (QD) hosted in 
a superconducting (S) substrate 21], In both SET and 
STM configurations such boson mode can be eventually 
related to external ac field. 


In what follows we introduce the Hamiltonian and dis¬ 
cuss the method for treating the bosonic mode. We next 
investigate the bosonic signatures in the QD spectrum 
and in the subgap Andreev conductance. For clarity, we 
focus on the limit -C ujq whereas the second coupling 
r$ can be arbitrary. In the last section we address the 
correlation effects. 
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II. MICROSCOPIC MODEL 


the molecular Hamiltonian ([2]) is transformed to 


For microscopic description of the tunneling scheme 
shown in Fig. [I] we use the Anderson impurity model 

H = H N + H s + H mol + H T . (1) 

Hn(S) refers to the normal (superconducting) lead, H mo i 
describes the molecular quantum dot (i.e. the local¬ 
ized electrons coupled with the boson mode) and Ht 
is a hybridization between the QD and itinerant elec¬ 
trons. We treat the normal electrode as a free Fermi 
gas H n = Dk.o- £kNCk aN Ck<jN and describe the other 
superconducting lead by the BCS Hamiltonian Hg = 
X)k,cr ^kSCk^CkcrS —A X^k (^kts^-k^s + ^-H^ki-s) ■ The 
annihilation (creation) operators correspond to mo¬ 
bile (3 = N,S electrons with spin a =t,4- an d energy 
£k /3 = £k /3 ^ 9p measured with respect to the chemical 
potential p,p. Nonequlibrium conditions can be driven 
by the bias V = /il — 9r and/or temperature difference 
Tj y T r . The induced currents depend qualitatively on 

the hybridization Ht = XX a p (^k/3 d^Ckcr /3 + H.c.^ and 
on parameters of the molecular quantum dot 

Hmol ^ E n\ a + U n dt n d i + ujq^ a + A ^ h dtJ (a 1, + a). 

G G 

The number operator h da = d' a d a counts the localized 
electrons with spin a, e is the QD energy level and U de¬ 
notes the Coulomb potential between opposite spin elec¬ 
trons. The boson field (described by fiW operators) is 
assumed as a monochromatic mode uio and its coupling 
with the QD electrons is denoted by A. 


H mol — E id a d a + Uh\n -|- + wofitfi, (3) 

G 

where the energy level is lowered by the polaronic shift 
e = e—A 2 /u>o and the effective potential U = U — 2A 2 /ujq. 

Boson operators are shifted d ' * = fiT) _ A d\d a 
whereas fermions are dressed with the polaronic cloud 

= d { J ] X m , X = e -(V^o)(a t -a)_ (4) 

Reservoirs Hp are invariant on the unitary transforma¬ 
tion © but the operator X appears in the hybridization 

term Ht- For simplicity we absorb it into the effective 
coupling constants Tp = 27tXX |Vk/3| 2 {X^X)6{uj — ^p) 
which can be defined for the wide band limit. 

The effective single particle excitation spectrum is 
given by the Green’s function 

Go(T\,T2) = —i (r T d a (Ti)dl(T2)^^ , (5) 

where T t denotes the time ordering operator. Since trace 
is invariant on the unitary transformations (...)^ = (...)^ 
it is convenient to compute the statistical averages with 
respect to H. In particular, © can be expressed as 

G< 7 (t 1 ,t 2 ) = -i (f T dcr(Ti)dl(T2)) - (f T X (ti) X* (r 2 )) - (6) 

because the fermionic and bosonic degrees of freedom 
are separated by the Lang-Firsov transformation. From 
a standard procedure d, |23[ one obtains 


III. MULTIPLE SUBGAP STATES 

There are three main obstacles in determining the ef¬ 
fective energy spectrum and the tunneling transmission 
of our system: i) the electron-boson coupling A, ii) the 
proximity induced on-dot pairing (due to A), and iii) the 
correlation effects caused by the Coulomb repulsion U. 
The most reliable way for studying them on equal footing 
would be possible within the numerical renormalization 
group [© approach, however such method encounters 
problems in estimating the Andreev transmission. To 
get some insight into the spectrum and transport prop¬ 
erties we start by neglecting the correlations and then 
(in the last section) treat them using the superconduct¬ 
ing atomic limit solution. 

Following [H-Qjl, we apply the unitary transforma¬ 
tion e s He~ s = H to decouple the electron from boson 
quasiparticles. With the Lang-Firsov generating opera¬ 
tor [ 23 ] 

s = — fi da (o^ - &) (2) 


(f T A(Ti)A(r 2 ) t ) = exp {-(A/w 0 ) 2 x (7) 
[(! - e" i “°( Tl - T2) )( 1 + N p ) + (1 - e“ o(Tl - T2) )A p ]} 


with the Bose-Einstein distribution N p = [e^ 0 — l] 1 . 
Fourier transform of the Green’s function © is found as 

G a (u>) = 0 ) e -(Vw^V"°) 2 (8) 


l 

xe'^o/ 2 /; 



Wo 


fy/N p (l + N p ) 


where /; denote the modified Bessel functions and 
9 a( T U T 2 ) = -i (T T d<T(ri)dt (r 2 )) a is the fermionic 

' ' H fer- 

part of ©. In the ground state © simplifies to 

hm G a (w) = J2g a (u}-luj 0 ) e~ 9 j 
1 

with the adiabatic parameter g = (A/wq) 2 . 


( 9 ) 
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Due to the proximity induced on-dot pairing the single 
particle Green’s function Gf(Ti,T 2 ) is mixed with the 
(anomalous) propagator 


F(ti,t 2 ) = —i ^T T d|(n)d|(T^)^ = (10) 

-*(t t ^(ti)4(t 2 )V (f T Xt( Tl )lt( T2 )\_ . 

' ' Hfer ' ' Hbos 


This important fact has been remarked in the previous 
considerations of dc Josephson current >1] and it also 
plays significant role for the Andreev spectroscopy (see 
the next section). The boson part of the anomalous prop¬ 
agator m takes the following form 


(f r X t (ri)X(r 2 ) t \ = exp {-(A/w 0 ) 2 x (11) 
[(1 + + jVp) + (1 + e iui °( Tl - T2) )N p \} . 

At zero temperature its Fourier transform simplifies to 

hm F(u>) = ^/(w-Zwo) e~ 9 (12) 


As regards the fermion part /(ti,t 2 ) = 

—i(T T 4( T i)4( r2 ))^ it couples to the Green’s 

function g(T\. r 2 ). Their Fourier components obey the 
Dyson equation 

g(uj) /(w) 

./*(-“>) -fl*(-w) 

- H° QD ( U ) - 


W — £ 0 

0 oj+s 


(13) 


where Y,° d is the selfenergy matrix of uncorrelated molec¬ 
ular dot and the second contribution is due to the 

effective Coulomb interaction U. In the wide-band limit 
the selfenergy Xq D (cc) can be expressed as 


^Qd( w ) 





(14) 


with 


7(w) 


VA 2 -lo 2 
i |u<| 

vT 2 -A 2 


for |w| < A, 
for |w| > A. 


(15) 


We investigated the effective spectral function p(w) = 
—7r _1 ImG(w + i0 + ) at zero temperature, focusing on the 
intermediate electron-boson coupling g ~ 1. Figures 00 
show the QD spectrum for 17 = 0, neglecting the correla¬ 
tion effects H c d rr . Influence of the Coulomb potential U 
is discussed in section V. 

Fig [2] illustrates evolution of the bosonic features with 
respect to the superconductor gap A. In the normal state 
(for A = 0) such lorentzian peaks are located at to = e + 
loj o (with integer l > 0) and their broadening is Idv + IV 
For finite A all peaks split into the lower and upper ones 



FIG. 2: (color online) Energy spectrum p(lo) of the uncor¬ 
related quantum dot (U = 0) obtained at T = 0 for i = 0, 
g = 1, u>o = 10rjv. The filled triangles at Iuj o ± A are only 
guide to eye. For increasing A the boson peaks split into the 
lower and upper states and their broadening shrinks to Fat. 


p(ca) 



FIG. 3: (color online) Spectrum of the uncorrelated quantum 
dot for ui o = 10F jv, g = 1, Ts = 4 Fjv, T = 0. The neighboring 
boson peaks are crossing at w = (| + l)uo for e — wo/2. 


due to the induced on-dot pairing. In the extreme limit 
A>Ts the selfenergy Sq D (uj) becomes static 

lim E° (u,) = -\( I s ) (16) 

A »r s QDK ’ 2 \ iT n ) K ’ 

therefo re the effecti ve quasiparticle energies evolve to 
loj o± sje 1 + (r s /2) 2 and their broadening shrinks to Tat. 

Focusing on such superconducting atomic limit USD 
we show in Fig. [I] the subgap bosonic peaks with respect 
to e. In the SET configuration the energy level e would 
be tunable by applying the gate voltage. In particular, 
these peaks may overlap with each other when e ~ wo/2 
as reported earlier in the Refs [H, Qjl, |20( . This effect 
can be deduced analytically from 

Icoo + Ve 2 + (r s /2) 2 = I'co o - Ve 2 + (r S /2) 2 . (17) 
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FIG. 4: (color online) Spectrum of the uncorrelated quantum 
dot for e = 0, g = 1, u>o = IOTat, T = 0, A 3> r.g. The 
bosonic features cross each other at ui = (^+Z)wo for Tg = too- 

The neighboring peaks (l' = 1 + 1) overlap when wo/2 = 

i Je * 1 + (T s /2) 2 . For small Tg such situation takes place 
at e ~ o- Other crossings would be eventually possible 
for the higher-order multiplications of ojq/2. 

Figured] displays the subgap spectrum p(ui) as a func¬ 
tion of the coupling Ts- From mu we conclude that for 

i = 0 the bosonic peaks overlap at Tg = ujq. Energy of 
these crossing points is w = (| + 1)ujq. Here (for g = 1) 
we observe four such crossings, but for stronger electron- 
boson couplings a number of the in-gap states and their 
crossings would increase. 

IV. ANDREEV CONDUCTANCE 


FIG. 5: (color online) The differential Andreev conductance 
Ga(V) versus the source-drain voltage V and the QD level 
e (tunable by the gate voltage). Results are obtained for 
T = 0, g s# 1, rs/Fjv=6, cjo/Fjv = 10, U = 0 and A 3> Tg. 
Conductance is expressed in units of 4e 2 //i. 

Fig. [5] shows the Andreev conductance as a function 
of voltage V applied between the metallic and super¬ 
conducting electrodes. We notice the differential con¬ 
ductance enhancements whenever V coincides with the 
in-gap quasiparticle energies. Since T A (— ui) = T A {u>) we 
observe these maxima at ±^/e 2 + (Tg/2) 2 ± Iloq. They 
eventually overlap when o is satisfied. In particular, 
for Ts = and ujq = 10Tjy the nearest bosonic peaks 
overlap when i = 4Fjv. Figure [5] clearly shows that the 
resulting maxima appear at |eV| = loq( 1/2 + l). 


Under nonequilibrium conditions the charge current 
can be transmitted at small voltage |eV| < A via the An¬ 
dreev scattering, engaging the in-gap states. This anoma¬ 
lous transport channel occurs when electrons from the 
metallic lead are converted into the Cooper pairs (prop¬ 
agating in superconducting electrode) with the holes re¬ 
flected back to N electrode. The resulting current Ia(V) 
can be expressed by the Landauer-type formula [24] 

Ia(V) = Jdu T a (ui) [f FD (u)-eV)-f FD (u+eV)}(, 18) 

with the Fermi-Dirac function f FF >(ijj) = [e aj / fc s T _|_ i] 
and the Andreev transmittance [241 ] 

T A (u)=T% |F(w)| 2 . (19) 

Optimal conditions for this subgap transmittance occur 
when oj coincides with the subgap quasiparticle states. 
In our present case we thus expect a number of such 
enhancements due the bosonic features. Let’s remark 
that T a (—u>) = T a (u}) implies the Andreev conductance 
G a (V) = 01 to be an even function of the bias V. 


V. CORRELATION EFFECTS 

In various experimental realizations of the quantum 
dots (such as self-assembled InAs islands carbon 
nanotubes [26j,[23] or semiconducting nanowires (28l.[29|) 
attached to the superconducting leads the energy gap A 
was safely smaller than the repulsion potential U. For 
this reason, in the subgap Andreev spectroscopy the cor¬ 
relations hardly contributed any Coulomb blockade. In¬ 
stead of it, they can eventually induce the singlet-doublet 
quantum phase transition [35[ and/or the Kondo physics 
[34| . In this paper we consider the strongly asymmet¬ 
ric coupling Tat <C Tg and focus on the deep subgap 
regime Tn,s ^ A, therefore the Kondo-type effects [30l — 
|34[ would be rather negligible. 

Analysis of such singlet-doublet transition for the vi¬ 
brating quantum dot has been previously addressed [u| 
using the NRG technique. We revisit the same issue 
here, determining the differential Andreev conductance 
(unavailable for the NRG calculations [l5[ ), because this 
quantity could be of interest for experimentalists. For 
the sake of simplicity, we analyze the correlation effects 
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G a (V) [4e 2 /h] 



FIG. 6: (color online) The differential Andreev conductance 
Ga(V) versus the Coulomb potential U and bias V obtained 
in the superconducting atomic limit A > Fs for T = 0, 
Fs/Fjv = 20 in absence of the boson mode <7 = 0. The thick 
(red) line indicates the QPT at [/ = Ts. 


FIG. 7: (color online) The subgap Andreev conductance 
Ga(V) as a function of the Coulomb potential U and volt¬ 
age V obtained for g = 1, ojo/Tn = 10 and the same model 
parameters as in figure [H] 


in the superconducting atomic limit A > Tj. Hamilto¬ 
nian of the molecular quantum dot ([3]) can be addition¬ 
ally updated with the pairing terms (d\d\ + d^dij 
originating from the static off-diagonal parts of the self¬ 
energy matrix (HI). 

In absence of the boson field (i.e for A = 0) the exact 
solution of such problem has been discussed by a number 
of authors (e.g. see the references cited in [36}). The ef¬ 
fective quasiparticle energies are given by ±{7/2 ± E d , 
where Ed = \J(e + U/2) 2 + (Vs/ 2) 2 . In the realis¬ 
tic situations only two branches ± (U/2 — Ed) appear 
in the subgap regime, whereas the other high energy 
states ± (U/2 + Ed) overlap with a continuum beyond 
the gap. The quantum phase transition (QPT) from the 
singlet u |0) ± v Iff) to doublet \a) configuration occurs 
at U/2 = Ed [35| . In order to estimate quantitatively 
the Andreev conductance we use the off-diagonal Green’s 
function f(to) [35, .36}, restricting to its subgap part 

j, f ^ a uv a uv 

fsub{0j) ~ u + Z£--(%-E d ) ~ io+^ + (^-E d ) 

( 20 ) 


with the usual BCS coefficient uv = Ts/^Ed and 


the spectral weight a = 


exp 


\2k B T ) 


exp 


\k B T J 


/z, 


where Z = 2 exp (±exp (^) ±exp (^) . The 
missing part of spectral weight 1 — a belongs to the high- 
energy states (outside the gap). At zero temperature this 
subgap weight changes abruptly from a = 1 (in the sin¬ 
glet state when U/2 < Ed) to a = 0.5 (in the doublet 
state when U/2 > E d ). 

In figure [6] we plot the Andreev conductance obtained 
for the half-filled quantum dot e = —U/2 (QPT occurs 
then at U = Ts). We notice the subgap conductance en¬ 
hancements around |eV| = U/2 — E d - Yet, exactly at the 
QPT, both the singlet and doublet contributions cancel 
each other. Formally, this is due to the odd (asymmetric) 
structure of the Green’s function (BUI) . 


The superconducting atomic limit solution can be gen¬ 
eralized onto g / 0 case in a straightforward way. The 
unitary transformation © implies £—>■£,[/—>■{/ and 
following the steps CEDED we can determine the off- 
diagonal Green’s function. At zero temperature, we find 


E sub((v) 


^ \ e” 9 (-«?)' /l\ 

auv } < - -f-z -r- 

^{uj+^-^-Ed) 

e ~ 9 (-gY /» 1 

_|_ ^E — Ed'j -s Iojq j 


with s = sign(^ — E d ). 

Figure [7] shows the Andreev conductance obtained for 
the half-filled quantum dot using g = 1, ojo/Fn = 10, 
r s /rAr = 20, T = 0. The bosonic side-peaks give rise to 
additional subgap branches, similar to what has been re¬ 
ported for the spectral function [ltj. Right at the QPT, 
the zero-bias conductance again vanishes G(0) —> 0 and 
we observe only the higher order maxima at |eF| = Iujo 
(with l > 1). Away from the QPT, the Andreev conduc¬ 
tance shows the usual maxima at |eF| = \U/2 — E d \+loJo 
whose spectral weights depend on U and l. 


VI. SUMMARY 

We have investigated the subgap spectrum and trans¬ 
port properties of the quantum dot coupled between the 
metallic and superconducting electrodes in presence of 
the external boson mode ojq. We have found that the in¬ 
duced Andreev states eventually cross each upon varying 
the gate potential (through e) or due to the correlations 
(via quantum phase transition from the singlet to doublet 
configurations). We have explored their signatures in the 
measurable charge transport. The tunneling conductance 
of such multilevel ’molecule’ shows a series of character¬ 
istic enhancements, dependent on: the gate voltage with 
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frequency ujq/2 (which can be deduced from Eqn. fTTl) . 
the bias V applied between external leads (Fig. [5]), and 
the correlations (Fig. [TJ) . External boson reservoir can 
thus substantially affect the anomalous Andreev current 
and it can be probed experimentally using the low-energy 
vibrational modes or the slowly-varying ac fields jlj. 
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